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ABSTRACT 

The origins of the bulge and disc components of galaxies are of primary importance to 
understanding galaxy formation. Here bulge-disc decomposition is performed simulta- 
neously in B- and _R-bands for 922 bright galaxies in 8 nearby (z < 0.06) clusters with 
deep redshift coverage using photometry from the NOAO Fundamental Plane Survey. 
The total galaxy colours follow a universal colour- magnitude relation (CMR). 

The discs of L* galaxies are 0.24 magnitudes bluer in B — R than bulges. Bulges 
have a significant CMR slope while the CMR slope of discs is flat. 

Thus the slope of the CMR of the total light is driven primarily (60%) by the 
bulge-CMR, and to a lesser extent (40%) by the change in the bulge-to-total ratio as a 
function of magnitude. The colours of the bulge and disc components do not depend on 
the bulge-to-total ratio, for galaxies with bulge-to-total ratios greater than 0.2. While 
the colours of the bulge components do not depend significantly on environment, the 
median colours of discs vary significantly, with discs in the cluster centre redder by 
0.10 magnitudes than those at the virial radius. Thus while star formation in bulges 
appears to be regulated primarily by mass-dependent, and hence presumably internal, 
processes, that of discs is affected by the cluster environment. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters 
- galaxies: bulges - galaxies: clusters - galaxies: photometry 



1 INTRODUCTION 

It is well-known that, at low redshifts, the centres 
of galaxy clusters are dominated by red, early-type 
galaxies. The morpholosy -density relation dates back to 
iHubble &: HumasorJ lll93lh . and was more accu r ately quan- 
tified bv lDressleJ (|l980i rand |Postman fc Gelled (| 19841 ). The 
early-type gala:xies (ellipticals and SOs) follo w tight scaling 
relations, such as the Fundamental Plane (jDressler et al.l 
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119871 : iDiorgovski &: David Il987 ) and the colour- magnitude 



relation (hereafte r CMR, ISandage fc Visvanathai 



lagnitude 
ajj 1 19781 : 



iBower et al.lll992l ) . Indeed it is the existence of a tight CMR 
for red gala^xies that allows one to isolate a "red-sequence" 
population. 

The dependence of the red-sequence CM R on 

environme nt is weak JSandage fc Visvanathan Il978ll . 
but real ([Abraham et aU 19961: Ivan Dokkum et al.l 19981: 



Balogh et al 



2000l: 
et al.l 



Carter et al. 2002: Lopez-Cruz et al 



2004: Bal ogh et al.l [20041 : iHogg et al.l |2004| : ISarazza et ~ 
2009.1 although in detail the strength of the effect depends 
on the morphological selection and on the magnitude range 
of the sample. These studies are based on global galaxy 
colours, from which it is impossible to disentangle whether 
such trends might be due to, for example, the colour of the 
bulge, the colour of the disc, or the varying relative impor- 
tance of bulge and disc components. 

Recent spectroscopic studies of the centres of red- 
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sequence giant galaxies have shown that, at a given mass, 
those in the cores of rich clusters are slightly older (by 
a few Gyr at most) than their counterparts in lower- 
density e nvironments, suc h as the outskirts of clusters 
Smith et al. 20061. g roups jProctor et al.l[2004l ') or the field 



Thomas et al., ,20051 : iBernardi et al.ll2006i l. There is little 



evidence of an environmental depe ndence of metallicit y (at 
fixed mass), although the results of ISmith et al.l ((2006) sug- 
gest that giant (~ L* ) red galaxies on the outskirts of clus- 
ters are ~ 15% younger and less enhanced in a-elements 
than those in the cluster cores. For fainter dwarf galaxies 
(Mr > —17 . 5), the enviro nmental dependence is stronger: 
ISmith et akl l|2008bl . [2009bl ') have shown that the stellar ages 
of dwarf red-sequence galax;ies in Coma vary by a factor ~ 3 
from the core to the virial radius. 

Spiral galaxies in clusters, when compared to their coun- 
terparts in the field, are re dder (Holmbcrg 1958; Kcnnicutt, 
119831: iBamford et al. 2009) and "anaemic" in appearance 
I van den Berghl 1 197& : .Mcintosh et al.l 2004), and poorer 
in ne ut ral hydrogen (|Davies fc Lewis! 19731: Sullivan et al.l 



19811: iH avncs et al' 'l984'; 'Giov anelli fc HavnesI 1 19851 : 



Cavatte et a l. 1990 ). Several mechanisms have been pro- 



posed for the dependence of spiral galaxy properties 
on environment: ram -pressure stripping of the cold gas 
Gunn fc GottI 1 19721 ). the removal of the hot gas halo 



Larson et al.' 1980), sometimes known as "strangulation" 



.,Balogh fc Morris 2000) and tidal st ripping by the cluster 
potential JMerrittj 1984 MamorJll987l) or by encounters with 



other galaxies (JGallaeher fc Ostrikeii 1 19721 : iMoore et al.l 
1 19981 ). The former mechanisms remove the gas that is the 
fuel for star formation, and so the stellar disc will dim and 
redden as its stars age. This will obviously affect the colour 
of the disc as well as the B/T ratio. In the tidal stripping 
scenario, outer parts of the stellar and gaseous discs are 
removed. A recent review of the effects o f such processes on 
late-ty pe galaxies in clusters is given bv lBoselli fc Gavazzil 
(l2006l) . 

From observations at higher redshifts, it is known that, 
in the past, clusters contained a higher proportion of spirals 
compared to SOs (Dressier et al. 1997; Fasano et al. 2000) 
than present-day clusters. Ivan Dokkum et al.l ([1993) have 
shown that SOs are bluer and have a larger scatter in colour 
than ellipticals in a z ~ 0.3 cluster, suggesting that some of 
them may have recently arrived on the red sequence. 

Clearly it is of interest to study the colours of bulge 
and disc components separately, since the separate com- 
ponents likely have different stellar ages and metallicities, 
and will be affected differently by environmental processes. 
Most previous studies have focused on the global colours 
of galaxies (sometimes subdivided by morphological class). 
Since most galaxies have both a bulge and a disc component, 
this global approach makes it difficult to disentangle the ef- 
fects of the two components. There have been few previous 
studies in which the separate colours of bulge and disc com- 
onen ts have been analyzed. In t he field. iBalcells fc Peletieij 






19941): iTerndrupet all (| 19941 ): iPeletier fc Balcellsl l|l99^ 



I Mac Arthur et al.l (200J), among others, have studied colours 



of disc and bulge components in spirals and SOs. In clus- 
ters, two-dimensional decompositions have been performed 
JCaon et al.l Il 990l: 'P'Onofric j [20011: ICutierrez et all |2004| : 
IChristlein fc Z abludofi 200l: Ide Jong et al.|[200i ). but the 



simultaneously is IKoo et al.l ((20051), who studied the z ~ 
0.83 cluster MS 1054-03. 

The goals of this paper are twofold. The first goal is to 
determine the origin of the slope of the red-sequence CMR 
by decomposing cluster galaxies into their constituent bulge 
and disc components, and studying the colours of each com- 
ponent separately. The second goal is to determine how en- 
vironment affects each of these components separately, as 
well as how it affects the global morphology, quantified by 
the bulge-to-total ratio (hereafter B/T). Other structural 
parameters (bulge effective radii and surface brightnesses, 
Sersic indices and disc central surface brightness and scale 
lengths) will be examined in future work. 

The outline of the paper is as follows. Section [2] de- 
scribes the sample of cluster galaxies, and details of the 
photometric and redshift data. Section [3] describes the mor- 
phologies and bulge-to-total ratios. In Section |4l the total- 
light CMR, as well as the CMRs for bulge and disc compo- 
nents is measured and its physical interpretation is discussed 
in Section [5] Section [6] quantifies the environmental depen- 
dence of the CMRs and the astrophysical implications are 
discussed in Section [7] . Throughout we adopt a Hubble pa- 
rameter, h — _ffo/(100 km s~^ Mpc^^) = 1 and a fim = 0.3, 
JIa =0.7 cosmology. 



2 DATA 

2.1 Cluster Sample 

As noted above, the goals of this paper are to determine the 
physical origin of the colour-magnitude relation and its envi- 
ronmental dependence within clusters. In particular, we are 
interested not only in the dominant red-sequence galaxies 
but also the contribution of blue galaxies. Photometry for 
the ~ 90 most luminous nearby clusters has been obtained 
as part of the NOAO Fundamental Plane Survey (here- 
after NFPS). The spectroscopic component of the NFPS 
(ISmith et al.ll2004l ) was, however, designed to follow-up only 
the red galaxies. In this paper, therefore, we focus on NFPS 
clusters with complete photometry for which nearly com- 
plete deep redshift samples are available from the litera- 
ture. Eight NFPS c luste rs previously studied by EN ACS 
(|Katgert et al.lll998l ). bv IChristlein fc Zabludo^ (|2003l ). by 
the C AIRNS project jRines et al.l |2003| ) and in the SDSS 
DR2 (jAbazaiian et al.l |200J) were selected, and are listed 
in Table [l] The typical cluster in our sample has a red- 
shift of 14 800 km s~^, a velocity dispersion, CTcI of 860 
km s~^, and a virial (7-200) radius of 1.46 fc~^Mpc, wher e 
7-200 = \/3(aei/1000km s-^)h~^M^c jCarlberg et al.lll997l ). 
The data extend to the virial radius, r2oo, and are complete 
within r2oo/2. Galaxies are assumed to be cluster members 
if their velocities lie within ±3(7ci of the cluster mean. De- 
tails of the galaxy completeness as a function of magnitude 
and colour will be discussed in Section [2.41 



only paper to perform this decomposition in multiple bands 



2.2 Field Sample 

In Section|6l we will also compare properties of cluster galax- 
ies to those in the field. In the foreground and background 
of the eight clusters, there are 35 field galaxies (after exclu- 
sion of two background clusters behind A0085 and A0151). 
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The median _R— band magnitude of the field sample is —20.3, 
very similar to that of the cluster sample. These field data 
are analyzed in exactly the same way as the cluster sample 
described above. 



2.3 Photometric Data 

The photometric data used were obtained as part of the 
NFPS (Smi th ct al. 2004 ), from the Kitt Peak National Ob- 
servatory (KPNO) 0.9 meter and the Cerro Tololo Inter- 
American Observatory (CTIO) 4.0 meter telescopes. The 
fields of view were 59' x 59' and 37' x 37' for KPNO and 
CTIO, respectively, and were centreed on the nominal X-ray 
centroid of the cluster. Exposure times were typically 400 
(R) and 600 (B) seconds at KPNO, and 60 (R) and 150 (B) 
seconds at CTIO. The seeing point spread function (PSF) 
full- width at half maximum typically was ~ 1.5 arcseconds 
for KPNO and ~ 1.1 arcseconds for CTIO. 

Data were reduced with the IRAF mscred package, and 
were cal i brate d to the Johnson-Cousins BR system using 
iLandoltl (|l992r ) standard stars. Object s and masks were gen- 
erated with SExtractor l|Bertin fc Arnouts,, 199(1 ). 



2.4 Corrections and Completeness 

Redshifts of galaxies within our images we re gathered from 
the NFPS spectroscopic survey (jSmith et a l. 2004) and from 
the NASA Extragalactic Database (NED). Galaxies were 
assigned cluster membership using an iterative Sa-clipping 
technique. Apparent magnitudes were converted into abso- 
lute magnitudes using cluster redshift in the Cosmic Mi- 
crowave Background (CMB) frame. 

Galactic exti nction corrections we re applied using the 
E{B-V) map of lSchlegel et al.l (|l998l . hereafter SFD), and 
converted to extinction using 4.325 E{B — V) for the B band 
and 2.634: E{B-V) for the i?-band (SFD). K-corrections 
were applied to the bulge and disc components separately, 
by interpola ting as a function o f morphological type from 
the tables in lFrei fc Gumil l| 19941 ). 

Fig. [T] shows that in general, the redshift data are com- 
plete at the bright end {Mr < —20), but that, as ex- 
pected, the completeness decreases at fainter magnitudes. 
To measure the completeness, a Schechter luminosity func- 
tion was fit to each cluster with ga laxies with Mr < —21 
fixing Mr = -21.14 and a = -1.21 (|Christlein fc Zabludofj 
I2OO3I ). but leaving the normalization, 0*, free. The fit is 
shown by the dotted line in Fig. [T] 

Based on Fig. [U we adopt a Galactic-extinction- 
corrected absolute magnitude limit of Mr < —19.5 for our 
primary sample. For the sample as a whole, the average com- 
pleteness brighter this limiting magnitude is 87%, and the 
completeness in the faintest 0.5 mag bin is 66%. To cor- 
rect for incompleteness, galaxies with redshifts are assigned 
a weight that is the inverse of the probability that a clus- 
ter member of that magnitude has a measured redshift. Note 
that the clusters A0085, A0119 and A0496 have redshift data 
that are considerably deeper. These three clusters are used 
to define a deeper Mr < — 18 sample, which will be used 
later to assess correlations at fainter magnitudes. 




Figure 1. i?-band magnitude histograms for individual clusters. 
The thin black histogram is all galaxies with photometry, whereas 
the thick red histogram shows those with published redshifts. The 
luminosity functions (normalized at the bright end by assuming 
a Schechter function) are shown by the blue dotted line. In each 
panel, the vertical dashed Une indicates the absolute magnitude 
limit adopted for that cluster. 



2.5 GIM2D Morphological Parameters 

Bulge and disc deco mposition is based on the GIM2D 
morphology package (jSimard et al.l I2OO2I ) which analyses 
postage-stamp images along with bad-pixel masks, PSF im- 
ages, and error maps and fits a PSF-convolved model galaxy 
to the two-dimensional images. The GIM2D model consists 
of a infinitesimally-thin exponential disc and an axisymmet- 
ric Sersic bulge. The fit yields bulge and disc magnitudes, 
radii, position angles, and inclinations (for discs), or ellip- 
ticities and Sersic index n (for the bulge components). 

Note that we could have chosen to fit the B and _R-band 
images independently, but in that case, it would have been 
possible to obtain different values in the two bands for pa- 
rameters such as the radii and position angles of given com- 
ponent. In particular, because of, for example, systematic 
errors in the PSF, there can be a trade-off between the flux 
in the bulge and the flux in the disc components. If such sys- 
tematic errors were different in the two bands, the colours of 
the bulge and disc components would not be robust. Instead, 
we used GIM2D in the two-band mode, in which simultane- 
ous fits to the two images are made for the radii, position 
angles, bulge Sersic index n, bulge ellipticity and disc incli- 
nation, while bulge and disc fl uxes are free in each filte r over 
all pixels in the SExtractor (JBertin fc Arnoutslll996l ) seg- 
mentation mask. While errors of the type described above 
can still affect these fits, in general, one would expect them 
to be yield more robust colours for the individual compo- 
nents. 

A limitation of two-band mode is that the bulge ef- 
fective radii are assumed to be the same in both B and 
R bands, and likewise for the disc scale lengths. In reality, 
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Table 1. Cluster Sample. N^ is the number of galaxies with Mji < —19.5. Phot is the source of photometry: K=KPNO; C=CTIO. 



Cluster 


R.A. 


Decl. 


<cz> 


N, 


A^lim 


O-cl 


'"200 


Phot 




(J2000) 


(J2000) 


km s~^ 






km s~^ 


ft-^Mpc 




A0085 


00 41 50.4 


-09 18 11 


16392 


144 


-18.0 


905 


1.57 


K 


A0119 


00 56 16.1 


-01 15 18 


12958 


80 


-19.5 


747 


1.29 


K 


A0151 


01 08 51.8 


-15 25 12 


15769 


51 


-19.5 


838 


1.45 


K,C 


A3128 


03 30 23.5 


-52 32 24 


17934 


107 


-19.5 


849 


1.47 


C 


A3158 


03 42 56.6 


-53 38 02 


17542 


88 


-19.5 


954 


1.65 


C 


A0496 


04 33 37.7 


-13 15 43 


9717 


127 


-18.0 


722 


1.25 


K 


A0576 


07 21 24.2 


+55 47 02 


11527 


60 


-19.5 


871 


1.51 


K 


A3667 


20 13 33.4 


-56 50 35 


16585 


93 


-19.5 


1000 


1.73 


C 



disc scale lengths are likely to differ, although for the SOs 
and early-type spirals studied here the diff erence is quite 
small: |Noordcrrnccr fc van der Hulstl ([2003) find that the 
disc scale lengths in a sample of SO-Sab galaxies are larger by 
3.2 ± 2.5% in the B-band than in the 7?-band. We have sim- 
ulated the impact that this GIM2D fitting constraint may 
have on disc colours by fitting model galaxies with disc scale 
length assumed to be different from the true one. We find 
that for a 3% error in scale length, the colour is biased at a 
level of 0.01 mags. Thus this is a small effect. 

2.6 Data Comparison 

Before proceeding with the analysis of the data, a check was 
performed to ensure internal consistency of the photomet- 
ric data and GIM2D decompositions. Of the clusters in the 
sample, only A0151 was observed both at KPNO and at 
CTIO. Figures [2] shows the comparison of the A0151 data 
obtained from the different observatories, in the sense A = 
KPNO - CTIO. In the magnitude comparisons (top panels 
of Figure [2)|, the median offset is 0.007 magnitudes in the 
7?-band and 0.025 magnitudes in the B-band. 

The scatter is larger than the individual measurement 
errors quoted by GIM2D. This suggests that GIM2D er- 
rors are underestimated, which is expected since these re- 
flect only the error due to photon noise, and not e.g. 
seeing mismatches, sky subtraction etc. Therefore, to the 
quoted GIM2D magnitude errors, we add in quadrature a 
magnitude-independent extra "systematic" error. Based on 
the scatter in the top panel of Fig. (2] we find that this extra 
error is 0.05 magnitudes. 

The B — R colour comparison is shown in the second 
panel of Figure [S] The median offset between KPNO and 
CTIO is 0.02 magnitudes, which is acceptable for the pur- 
poses of this study. The measured scatter suggests that in- 
dividual colour measurements are accurate to 0.03 magni- 
tudes. This is smaller than the magnitude error in either 
band separately, and is due to the fact that the fit is made si- 
multaneously in both bands. From similar comparisons (see 
bottom panels of Fig. [2]), errors in a single measurement of 
the bulge and disc colours are estimated to be 0.07 and 0.11, 
respectively. 

Figure[3]shows the differences in the bulge-to-total light 
ratios. The observed scatter in the R-band B/T values is 
larger than the GIM2D errors. Therefore, we add in quadra- 
ture a systematic error of 0.07 per individual B/T measure- 
ment. Given the random (noise) errors, the total measure- 
ment error in B/T is then typically ~ 0.08. 



0.4 


r 




' Median: 
Scatter: 


d.ooT 

0.064 


' ' 1 




1 1 1 1 1 




--_ 


O.S 



- 




• 


" • • 


-•r.-'l 


.• 


• • 




- 






* • 






• 


• : •• V 




• 


0.2 


r 






. 










- 


0.4 


r 






1 1 1 


1 1 1 








- 


0.4 


r 




' Median: 
Scatter; 


(J.O34' ' 
0.045 






1 1 1 1 




- 


o.a 




r 




• 


-..- 'v.- 


•I'. .. 


;■ 




• 


* ~ 






I 




* • 








- 


0.3 


r 
















- 


0.4 


r 
















-. 


4 






'Median:' 


(iM-i ' 






1 




■ 








Scatter: 


0.091 . 






• 




: 


0.2 



- 




, 


• 


• * '. 


. 


• °' 


■ 


- 






"^ 






. 


''*.-• 




: 


0.2 


r 












• 




- 


0.4 


r 






1 1 1 


1 1 1 




• 




H 








'Median:' 


-A.oii ' 






1 




- 








Scatter: 


0.147 




, 


• 




; 


0.2 



0.3 


- 


. 




• 


••>" 




• 


° 


~ 


- 




• 


•' 


: • 


■ 




*'° 




0.4 


r 






1 


1 




1 




- 



Average Magnitude 



Figure 2. Comparison of the total magnitudes and colours as 
measured by GIM2D from independent runs. A is in the sense 
KPNO — CTIO. The upper panel compare the H-band magni- 
tudes, the second panel shows the global B — R colour compari- 
son, and the bottom two panels compare bulge and disc colours. 
The median offsets are small. 



3 MORPHOLOGIES AND BULGE-TO-TOTAL 
RATIOS 

3.1 Visual Morphological Classification 

For a subset of the sample described above, postage stamp 
images, pixel masks and GIM2D residual images for all 
galaxies were inspected visually by one of us (GAW) and 
a visual morphological classification was assigned to each 
galaxy. Furthermore, cases where the automated pixel mask- 
ing had failed (e.g. close galaxy pairs or missed bad columns, 
etc) were flagged and rejected from the sample. 

Fig.|4]compares the R-band bulge-to-total (B/T) ratios 
from GIM2D with the visual classiflcation. The visual classi- 
flcations were labeled according to the corresponding T-type 
and plotted against B/T. The trend of decreasing B/T with 
increasing T-type, and the dispersion in B/T (~ 0.17) at a 
given type, ar e in genera l ly in good agreement with pre - 
vious studies IjKentl 1 19851 : ISimien fc de VaucouleursI 119861 ) , 



Colours of Bulges and Discs within Galaxy Clusters and the Signature of Disc Fading on Infall 5 



c 

IB 
I 



CQ 
< 



Slope: 0.077 +/- 0.153 
Scatter: 0.137 



Hu 



0.2 0.4 0.6 o.a 1 

Average B/T 

Figure 3. Difference between the ij-band B/T measured from 
CTIO and from KPNO as a function of average ij-band B/T. 
Tlie plotted GIM2D errors underestimate the differences, which 
presumably arise from other systematic sources of error such as 
modelling the PSF. 



although the median B/T ~ 0.3 for Sc galaxies is some- 
what higher than expected. Since the measurement error in 
B/T is typically 0.08 (Section 12. 6p . the intrinsic scatter in 
B/T of a given T-type is ~ 0.15. 



3.2 Tests of the GIM2D decomposition 

GIM2D fits a Sersic bulge plus disc model to the B— and 
7?— band data simultaneously. While we argue above that 
the constraints of these fits lead to more robust colours of 
the bulge and disc components, it's not necessarily true that 
these components — needed to reproduce the 2-D light dis- 
tribution — correspond to real bulge and disc components. 
For example, a pure bulge component with an isophote twist 
might be fitted with a (untwisted) bulge and a "disc" which 
accommodates the twist. 

One way to check the reality of the bulge and disc com- 
ponents is by plotting their ellipticity distributions. If the 
"discs" are really related to the bulge components, they 
would be expected to follow the ellipticity distribution of 
bulges, i.e. there should be few with e = 1 — b/a > 0.5. 
Fig. [5] shows the ellipticity distribution for bulges and disc 
components. If the discs are infinitely thin, transparent and 
seen at random inclinations, then they would have a uni- 
form distribution in e. The discs in Fig. [S] are very close to 
this scenario, although there is a deficit of very high e > 0.9 
discs. We attribute this to the finite vertical scale-height of 
discs: even seen edge-on, such discs will have b /a > 0. Note 
that t his result is different from that found bv lDriver et al.l 
l|2007h for discs in the field. In addition to the deficit of very 
high £ discs, they also find a deficit of moderately-inclined 
discs: for example, in the range 0.5 < e < 0.7 (corresponding 
to 60° <C J <; 73) their counts are reduced by ~ 40% (cf. their 
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Figu re 4. Visual morphologica l classifications encoded by T- 
type llde Vaucouleurs et al.lll99J) versus R-band B/T. The filled 
circles show the mean B/T at each T-type, and the long errorbars 
indicate the dispersion in B/T at each type. The short errorbars 
are the error in the mean. Small dots are individual data points, 
with the T slightly randomized for graphical purposes. The num- 
ber of galaxies in each bin is also indicated. 



Figure 5), whereas ours are consistent with being constant. 
They conclude that dust obscuration dims these moderately- 
inclined galaxies resulting in their under-representation in a 
magnitude-limited survey. In contrast, our uncorrected disc 
e distribution looks very similar to their distribution after 
correction for dust extinction. The difference between the 
field and cluster disc ellipticities suggests that cluster discs 
may be depleted in dust, and hence considerably more trans- 
parent. In contrast to the flat distribution of disc ellipticities, 
the bulge components are clustered between 0.1 < e < 0.5, 
as expected ba sed on the distribution of E-I-SO ellipticities 
in clusters (e.g. Ijorgensen fc Franxlll994h . 



3.3 Dependence of B/T on Magnitude and 
Environment 

Before discussing the colours of bulge and disc components, 
it is interesting to study the dependence of B/T on magni- 
tude and environment. Throughout this paper, unless explic- 
itly stated otherwise, B/T refers to the 7?-band value. Fig. [6] 
shows the weak, but statistically significant, dependence of 
the bulge-to-total ratio on magnitude, with more luminous 
galaxies being more bulge-dominated, mirroring the trend 
of morphological type with luminosity. 

The morphology-radius relation, where morphology is 
quantified by B/T and the projected radius, rp is scaled 
by the virial radius, r2oo, is shown in Fig. [T] The variation 
in B/T is rather weak: only in the very innermost regions 
(rp <^ 0.2 r 200 ) is there evidence for a difference in the me- 
dian B/T. However, a Spearmank rank correlation test does 
suggest a correlation at the 98.7% confidence level. More- 
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Figure 5. Distribution of cUipticity for disc (dashed blue) and 
bulge (solid red) components. Note that there is not an excess of 
disc components with ellipticities typical of bulges. The two dis- 
tributions are very close to what one expects a priori on physical 
grounds indicating that the disc and bulge components are not 
artifacts of the fitting procedure. 
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Figure 6. B/T in the ij-band as a function of magnitude for the 
Mjj < —19.5 sample. The square symbols with error bars show 
the median B/T in a bins of magnitude. The solid line is a fit to 
the medians. 
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Figure 7. B/T in the ij-band as a function of scaled projected 
radius for the Mji < —19.5 sample. The square symbols with 
error bars show the median B/T in a bin of rp/r2oo- The line is 
a fit to the medians. 



over, fitting B/T as a function of log(rp/r2oo) yields a signif- 
icant slope: — 0.108 ±0.034J. Thus while there is a clear dif- 
ference between morphologies in clusters and in the field, the 
morphological segregation as defined by B/T within clusters 
appears to be limited to the innermost regions. 

Finally, it is interesting to compare the cluster re- 
sults with those in the field. The distribution of B/T (for 
Tp < 0.5r-2oo) is shown in Fig. [S] For comparis on, we also 
show the B— band B/T distribution in the field (|Allen et al.l 
l2006f ) for the same absolute magnitude limit as applied to 
the NFPS clusters. Compared to the field, clusters contain 
an excess of galaxies with SO morphology (0.4 < B/T < 0.6), 
and a clear deficit of spirals, particularly those of the latest 
types, or lowest bulge fractions {B/T < 0.1). The distribu - 
tion of B/T in the group environment (jMcGee et al.l 120081 ) 
is intermediate between that in the field and in clusters. 



4 THE RED-SEQUENCE 

COLOUR-MAGNITUDE RELATION 
REVISITED 

In this section we examine the global-colour vs. magnitude 
relation for individual clusters as well as the universal re- 
lation for all clusters. We then decompose the global-CMR 
into the CMRs for bulges and discs separately. 



^ In Fig. 1141 below, we demonstrate that there is no significant 
magnitude-radius relation, hence the correlation found here is not 
an artifact of luminosity segregation coupled with the dependence 
of B/T on magnitude. 
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Figure 8. The solid red histogram shows the distribution of 
fl-band B/T for Mn < -19.5 galaxies in the NFPS clusters, the 
dashed blue line is the same but for the B— band B/T, and, for 
comparison, the dotted blue lin e shows the B— band B/T for the 
field sample 11 Allen et al.ll2006h . selected to the same magnitude 
limit as the cluster sample {Mn < —19.5). The most striking dif- 
ferences are the excess of late-type spiral {B/T < 0.1) systems in 
the field, and the excess of SO-like (0.4 < B/T < 0.6) morpholo- 
gies in the cluster. 



4.1 Robust Fits to Galaxies within r2oo/2 

We will show below (Section [6]) that the colours of some 
components depend on cluster-centric radius. Consequently, 
in this section, we restrict fits to galaxies with projected 
cluster-centric radii Vp < r2oo/2. Furthermore for fits to the 
red-sequence, we do not wish to be biased by very blue galax- 
ies, nor by outliers with anomalous colours. Therefore, in 
order to fit the CMR in a robust way, the data were first 
binned by magnitude, with an equal number of data points 
in each bin. For each bin the median was calculated, and a 
linear least squares trend was fit to these medians. 

To estimate the scatter of the data, the semi- 
interquartile range (SIQR) of the data is calculated, which 
is then converted into a the equivalent dispersion {a = 
1.5SIQR) of a Gaussian with the same SIQR. In order 
to compare the colour-magnitude relations with different 
slopes, a normalized colour is defined, (B — R), the value 
of the fit at M^ — —21.14. The errors on the normalized 
colour were calculated using 1000 bootstraps. 



4.2 Global Colour-Magnitude Relation 

The global colour-magnitude relation is shown in Figure |9l 
For galaxies with AIr < —19.5 the normalized colour at 
M^ is {B - R), = 1.485 ± 0.009 magnitudes, with a CMR 
slope of -0.031 ±0.006. For the fainter subsample of A0085, 
A0119 and A0496 combined, which extends to 7? < -18.0, 
the slope steepens slightly to — 0.042±0.00 7 with the normal- 
ized colour being 1.463 ±0.012. Note that iLopez-Cruz et al.l 



Figure 9. The colour-magnitude diagram for all 8 clusters com- 
bined. Upper panel: Filled circles indicate total-light colours of 
individual cluster galaxies with rp < 0.5r2oo- Median colours in 
bins of magnitude are shown by the open squares. The robust fit 
to the median colours is shown by the solid line is to galaxies with 
AIr < —19.5 and r'p < r2oo/2. Lower panel: Same median data 
as in the upper panel, but shown at a different scale to better see 
the slope. 



(!2004) find a CMR slope of -0.052 ±0.008 for Mr < -18.6, 
which is consistent with our result within the uncertainties. 
The "1-sigma" scatter in the colour-relation is 0.07 mag- 
nitudes, as estimated via the SIQR. Note that the mea- 
surement error in B — 7? is estimated to be 0.03, so this 
yields a intrinsic scatte r of 0.06. This is consisten t with the 
value (0.074) found bv iLopez-Cruz et all (|2004i ). although 
their limiting magnitude {Mr = —18.6) is somewhat fainter 
whereas their radial limit (0.2 r2oo) is smaller than ours 
(0.5r2oo). 



4.3 Fraction of Blue Galaxies 



The blue fraction was defined bv lButcher fc Qemleij (Il984l . 
hereafter BO) to be the fraction of galaxies brighter than 
A'lv = —20, lying within a third of the cluster radius and 
bluer than the red sequence by 0.2 mags in B — V. After 
correcting for cosmology, we find that their magnitude cor- 
responds to Mr — —19.1, which is only 0.4 mags fainter 
than our completeness limit. To match the BO definition 
as closely as possible, here we define a blue galaxy is as 
one with {B — R) < 1.25 after having normalized to M* 
by correcting for the CMR slope, i.e. 0.24 magnitudes bluer 
than the red sequence. Within 0.5 r2oo there are only 35 
blue galaxies of 595 in total in the 8 clusters, yielding a 
blue fraction of 5.9± 1.0%. For A0085, A0119 and A046, the 
clusters complete to Mr < —18, within 0.5r2oo, the blue 
fraction climbs from 12/183 (6.5 ± 1.9%) at Mr < -19.5 to 
27/259 = 10.4 ± 2.0% at Mr < -19. This suggests that the 
blue fraction for the whole sample to AIr < —19 is 10 ± 2%. 
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This blue fraction is considerably less than the blue fr a ctions 
(~25%) at the same radius quoted in lAndreon et alj (|2006|) 
for clusters at z ~ 0.35. 



4.4 Colour-Magnitude Relations for Individual 
Clusters 

The individual cluster colour-magnitude relations are shown 
in Figure [TOl The normalized colours (B — R), for each clus- 
ter are consistent with the global value. Only the slope in 
A0085 deviates by more than 2a from the global slope. Thus 
there is no evidence for deviation from the universal CMR 
in our sample. 



4.5 Bulge and Disc Colour-Magnitude Relations 

Having established that the total colours are uniform and 
consistent with previous results, we now turn to the main 
goal of this paper. 

Our two-band GIM2D fits yield the colours and CMRs 
of the bulge and disc components separately. However, 
galaxies with very low B/T do not have reliable bulge 
colours, so, in measuring the bulge colour-magnitude rela- 
tion, we have omitted galaxies with B/T < 0.2. This culled 
bulge colour-magnitude diagram is shown in the left panel of 
Figure [TT] Note that the horizontal axis is the total 7?— band 
magnitude. The slope of the bulge CMR, -0.037 ± 0.014, is 
very similar to the slope of total-light CMR, however the 
normalized colour, 1.581 ± 0.091, is 0.1 magnitudes redder 
than the total colour. Using the population synthesis models 
of Maraston, for an assumed metallicity [Z/Zq] = 0.3 dex, 
this colour corresponds to an SSP age of 12 Gyr. At fainter 
mags Mr < -18 the bulge CMR steepens to -0.067± 0.017. 

The measured scatter in the bulge CMR (0.093) is ac- 
tually somewhat larger than the total CMR scatter (0.07). 
This can be understood if one allows for additional error in 
the bulge-disc decomposition itself: some disc light is incor- 
rectly assigned to the bulge and vice versa. In Section [2l 
we estimated from repeat measurements that the measure- 
ment error on the bulge colour is ~ 0.07, and so the intrinsic 
scatter is likely to be much lower. 

If we fit bulge colour as a function of bulge magnitude, 
we obtain a flatter slope —0.025 ± 0.013, and larger scatter 
(0.098). However, in Section[S]we will show that bulge colour 
is most closely linked not to total magnitude but rather to 
central velocity dispersion. 

The disc colour-magnitude relation for galaxies with 
B/T < 0.8 is shown in the left panel of Figure [Til We 
find that discs have a flat CMR: the measured slope is 
O.OOOiO.OlO. The disc colour for an L. galaxy (1.353±0.016) 
is significantly bluer than that of the bulge component at the 
same total magnitude. Thus, in clusters, discs are bluer than 
bulges by 0.228±0.024 mags at L*. The bulge-disc colour dif- 
ference is larger than, but m a rgina lly consistent with, that 
found by iPeletier fc Balcellj (| 19961 ) for SO-Sbc galaxies in 
the field (0.045 ± 0.097). On the other hand, the colour dif- 
ference in clusters is smaller t han the value 0.301 zfc 0.014 
for a SO/a-Sdm field sample JMacArthur et all IJOoi ). At 
fainter magnitudes, the difference in colour is much less: at 
Mn = — 18, ~ 0.09 ± 0.06 mags. The colour-magnitude re- 
lations are summarized in Table (2] When disc colours are 
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Figure 12. Histograms of the colours of disc (solid lines, for 
galaxies with B/T < 0.8) and bulge (dotted lines, for galaxies 
with B/T > 0.2) colours, for subsamples with different B/T. 
Note there is little dependence of disc or bulge colour on B/T, 
except for discs in very late types (B/T < 0.2). 



Table 2. Galaxy colours and properties as a function of magni- 
tude. 



Property 



Mag Limit. 



slope 



value at L* 



Total Colour M < -19.5 

Bulge Colour M < -19.5 

Disc Colour M < -19.5 

B/T M < -19.5 



Total Colour 
Bulge Colour 
Disc Colour 
B/T 



M < -18 
M < -18 
M < -18 
M < -18 



-0.031±0.006 

-0.037±0.014 

O.OOOitO.OlO 

-0.069±0.020 

-0.042±0.007 

-0.067±0.017 

0.000±0.014 

-0.069±0.019 



1.485±0.009 
1.581±0.009 
1.353±0.016 
0.614±0.020 

1.463±0.010 
1.596±0.016 
1.307±0.029 
0.643±0.026 



regressed against disc magnitude the slope is slightly neg- 
ative, although the difference from zero is not statistically 
significant: —0.006 ±0.015, and the scatter is similar (0.120) 
to the scatter versus total magnitude.. 

It is important to investigate whether bulge and disc 
colours d epend o n morphology. For example, in the field, 
[MacArth ur et all ([2004) found that the colours of both 
bulges and discs are correlated with morphology in the sense 
that later morphological types host both blue bulges and 
blue discs. Figure [12] shows a histogram of the bulge and 
disc colours for different bins of bulge-to-total ratio. The 
bulge and disc colours show no significant dependence on 
morphology for B/T > 0.2. Discs in B/T < 0.2 galaxies are 
significantly bluer than the discs in higher B/T galaxies. 
However, there are few (9%) such low B/T galaxies in our 
sample. 
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Figure 10. The colour-magnitude diagrams for each of the clusters. The solid line indicates the least-squares fit for the data in the 
cluster, and the dashed line shows global colour-magnitude relation. The fits are to the Mjj < —19.5 sample. 



4.6 Deconstructing the Colour-Magnitude 

Relation: Contributions from Bulge and Disc 
colours and from the BuIge-to-Total ratio 



From the bulge and disc colour magnitude diagrams, it ap- 
pears that the global CMR slope is due primarily to the 
colour of the bulge component. A quantitative check of this 
conclusion can be performed by reconstructing the global 
slope using the bulge and disc CMRs and the dependence 
of B/T-raiio on magnitude. Specifically, the global colour- 



magnitude relation can be decomposed as 
d{B - R) _ dB/T 



AMr 



+ 



+ 



OMr 
d{B-R)B 

OMr 
d{B-R)D 



[{B-R}^-{B-R)^] 
(B/T) 
(1 - (B/T)) 



OMr 

For our Mr < —19.5 sample we measure an average 
{B/T) = 0.566 ± 0.04 and an average difference in bulge 
and disc colour is {B - R) g - {B - i?)^ = 0.237 ± 0.013. 
Inserting the measured slopes from Tableland evaluating 
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Figure 11. (Left Panel) Bulge colour vs. total H-band magnitude. All galaxies which are disc- dominated {B/T < 0.2) have been 
removed. Symbols are as in Fig. 1101 The solid line is the robust fit to bulge colours with Vp < r20o/2 and Mn < —19.5. The dashed 
line is the total-colour-magnitude relation for comparison. (Right) Disc colour vs. total R-band magnitude. All galaxies which are bulge 
dominated {B/T > 0.8) have been removed for this comparison. Symbols and curves are as in the left panel. 



each term yields -0.014 ± 0.004 for the first {B/T) term, 
-0.021 ± 0.008 for the second (bulge CMR) term and zero 
from the third (disc CMR) term. 

Thus the origin of the "tilt" or slope of the total-light- 
CMR is due mostly (60%) to the slope of bulge CMR, with 
the remainder of the total-light-CMR slope arising from the 
change in B/T as a function of magnitude, which mixes 
in increasing more bluer disc light at fainter magnitudes. 
These trends, however, do not give insight into the ages or 
metallicities of the stellar populations of our sample. We 
address this topic in Section [S] below. 



5 DISCUSSION: THE ORIGIN OF THE CMR 

We have shown that the tilt of the CMR is due primarily 
to the bulge CMR, and to a lesser extent, to the change 
in morphology along the red-sequence. The bulge-CMR it- 
self has been thought to be due to an age-luminosity rela- 
tion, a metallicity-luminosity relation, or both. While opti- 
cal colours alone cannot break the "age-metallicity" degen- 
eracy, in principle it is possible to do so by studying spec- 
troscopic absorption linestrengths. lSmith et al.l (|2009a|) and 
[Graves et al. (2009) have shown that the dominant "driver" 
of stellar populations is the central velocity dispersion, a, 
and that any dependence on, for example, stellar mass is 
secondary. 

It is therefore interesting to examine the colours as a 
function of velocity dispersion, and to model these colours 
as a function of stellar pop ulation parameter s. Using the 
NFFS velocity dispersions (jSmith et al.l l2004i ) for a sub- 
sample of emission-free, red galaxies with Mr < —19.5, 
we obtain the colour-cr relation (CSR) shown in Fig. 1131 



the slop e is 0.27 ± 0.04, consistent with previous determi- 
nations (JMatkovic fc Guzma"nll2005l ). The bulge-CSR has a 
statistically-significant slope of 0.21 ± 0.08, whereas that of 
the disc is statistically insignificant (0.09 ±0.09). We note in 
passing that the scatter in colour, for bulges and for discs, 
at fixed a is lower than the scatter at fixed magnitude (for 
the same red-selected subsample), suggesting that velocity 
dispersion and not, for example, stellar mass determines the 
colour of the disc as well as the bulge. 

The fiatness of the disc-CMR and the disc-CSR suggests 
that there is no significant variation in the stellar popula- 
tions of the disc component along the sequence. As shown 
in the Appendix, the disc colour is consistent with an SSP 
age of 4.5 Gyr and solar metallicity. 

It is possible to combine the disc colours, the GIM2D 
fits and nuclear linestrength data to derive ages and metallic- 
ities of the bulge component. These fits are described in de- 
tail in the Appendix, we summarize the results here: the fits 
require a variation in bulge age, in the sense that bulges in 
lower velocity dispersion galaxies are younger, and a weaker 
variation in bulge metallicity. 

An age variat ion along t he bulge sequence is consistent 
with the results of lSmith et a l. (2006) from the NFFS. They 
included B/T as an additional parameter in their fits of age, 
metallicity and a-enhancement, and found that bulges are 
somewhat older than discs, by, on average, a factor 1.5, but 
that the effect of B/T was weaker than the dependence on 
velocity dispersion. Thus, if discs are 4.5 Gyr old, one would 
expect bulges to have an age of ~ 7 Gyr. This is reasonably 
consistent with the fits described in the Appendix for which 
bulge ages ranging from 5 Gyr to 12 Gyr, as a function of 
velocity dispersion. The data are certainly consistent with 
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Figure 13. Colour versus log((T) for galaxies in the NFPS-II sample (red and emission free), with otherwise the same selection criteria 
as in Figures l9l and [TTI The three panels show the total colour, bulge colour and disc colour, from left to right. 



scenario in which the bulge CSR (and CMR) "tilt" is largely 
due to age, with metallicity being a secondary effect. 



6 ENVIRONMENTAL TRENDS WITHIN 
CLUSTERS 

It has long been known that galaxy properties are corre- 
lated with environment, the most wel l-established tr end be- 
ing the morphology-density relation (iDressleiJIlQSOf) . While 
the dependence of total colours on en vironment has been 
establis hed by a number of authors ([Balogh et al.l |200J; 
[Lopcz- Cruz et all |2004| : iHogg et al.1 120041 '). there has been 
little attention paid to the colours of bulges and discs seper- 
ately. In this section, we consider two indicators of envi- 
ronment: the projected cluster-centric radius, and the local 
density. The former is connected to the properties of the 
cluster as a whole, whereas the latter is more sensitive to 
local structure. 



6.1 Dependence on Cluster-centric Radius 

Before turning to the dependence of colours on measures 
of environment, we must first examine the possibility of lu- 
minosity segregation, because, if there were significant lu- 
minosity segregation in clusters, one might expect to see 
a colour-radius effect arising from the combination of lu- 
minosity segregation and the CMR. Fig. [14] shows magni- 
tude as a function of radius. The brightest galaxies, with 
Mji < —22 appear to prefer rp/r2oo < 0.6. For the remain- 
der of the galaxies, no evidence for luminosity segregation 
in our s ample. This result is in agreement with the recent 
studies jAdami et a l.'llQQS": iBiviano et"aLll2002l:|Pracy et all 
l2005l : lMercurio et a l. 2006; von der Linden et alBoiOl ). Fur- 
thermore, when if we subtract the CMR from the observed 
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Figure 14. i?— band absolute magnitude as a function of pro- 
jected cluster-centric radius. 



colours and look for the dependence of the colour residu- 
als on radius, we find trends of similar significance to those 
found below. Hence for simplicity, we simply fit colour as a 
function of cluster-centric radius. 

Figure [15] shows total, bulge and disc colours versus 
projected cluster-centric radius in units of the virial radius: 
rp/r2oo, Here, only the data with Mr < —19.5 are used. 
The left panel shows the total colour as a function of radius. 
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Figure 15. Total colour, bulge colour, and disc colour versus projected radius in units of r200 for the Mji < —19.5 sample. The upper 
panels show individual points and medians (open squares). The filled square denotes the median for the field sample, and is placed 
arbitrarily at 1.05 r/r2oo- 



The slope of the fit to median colours is -0.020±0.025. The 
bulge colours are shown in the middle panel. The slope of the 
bulge-colour-radius relation is not statistically significant. 
Finally, the disc colours are shown in the right-hand panel, 
where the radial slope is statistically significant at the 4. So- 
level: —0.103 ± 0.023. These fits, as well as those described 
below, are summarized in Table [S] 

Visual examination of Fig. [TS] may hint that much of 
the change in disc colours is due to a tail of blue discs in 
the outskirts of clusters. However, recall that we are fitting 
the median colour (not the mean) and this is much less sen- 
sitive to outliers. The disc colours for "core" (rp < 0.3r2oo) 
and "outskirt" (0.3 < rp/r2oo < 1) samples are shown in 
Fig. 1161 This confirms that, while such a blue tail does ex- 
ist, the peak of the distribution is also shifted to the blue 
in the "outskirts" bin. Moreover, this trend continues to the 
field: the peak of the colour distribution shifts blueward, 
and the fraction of very blue discs [B — R <1) increases. A 
Kolmogorov-Smirnov (K-S) test shows that the "core" bin 
is significantly different from both the "outskirts" bin and 
from the field, at 99.92% and 99.5% confidence levels respec- 
tively. Thus there is a clear difference between the cores of 
rich dusters and the field. 

We have also performed the same bulge and disc fits 
as a function of log(rp/r2oo) (see Table [31). This functional 
form gives greater leverage to the innermost regions. For 
this parametrization, the radial dependence of the total 
colour is highly significant: the slope is —0.092 ± 0.013 per 
decade in radius. Furthermore, with this parametrization, 
the bulge colour dependence is also statistically significant 
(-0.046±0.018), while the disc dependence (-0.065±0.015) 
remains significant. The difference between the results from 
the linear fit and the logarithmic fit arise from the "core" re- 
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Figure 16. Histograms of disc colour for two radial ranges rp < 
0.3r2oo (rsd solid) for rp > 0.3r2oo (green dashed) and in the field 
(blue dotted). The vertical marks indicate the median colours of 
each sample. 



gions rp <C 0.2r2oo, where bulge colours (and consequently 
the total colours) appear to be somewhat redder than at 
larger annuli. 

We have also attempted to extend this analysis to 
fainter magnitudes but found that, when the sample is re- 
stricted to only the three clusters complete to Mb. < —18, 
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Table 3. Galaxy colours and properties as a function of environ- 
ment. Slopes that differ from zero at greater than 3cr significance 
level are shown in bold type. 
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X 


slope 


Mr 


rp/r200 


0.063±0.194 


B/T 


rp/r 200 


0.004±0.061 


Total Colour 


i'p/r200 


-0.020±0.025 


Bulge Colour 


rp/r200 


0.007±0.038 


Disc Colour 


rp/r200 


0.103±0.023 


Mr 


log{rp/r20o) 


-0.182±0.113 


B/T 


log(rp/r2oo) 


-0.108±0.033 


Total Colour 


log(rp/r20o) 


-0.092±0.013 


Bulge Colour 


log(rp/r2oo) 


-0.046±0.018 


Disc Colour 


log(rp/r20o) 


-0.065±0.015 


Mr 


log(S5) 


0.054±0.064 


B/T 


log(E5) 


0.048±0.019 


Total Colour 


log{S5) 


0.032±0.008 


Bulge Colour 


log(S5) 


0.006±0.014 


Disc Colour 


log(E5) 


0.050±0.013 



7 DISCUSSION 

7.1 Bulges 

The bulge colours have show no dependence on environment 
from the cores of rich clusters to the field. At face value, this 
is compatible with the results fro m morpho logically-selected 
samples. For example, van Dokk um et al.l (^1998) found no 
clustercentric variation in their E sample (although they did 
find a trend for their SO sample) . iHogg et al.l (|200J) found a 
very weak dependence of the peak of the colour distribution 
of Sersic n > 2 gala xies with density. Finally, in a cluster 
at 2 = 0.83. IKoo et a l. (2005) also found that bulge colours 
were insensitive to cluster-centric radius. 

However, we note that it is possible for the stellar 
populations of bulges to change in such a way as to pre- 
serve B — R colour. For this to occur, changes in metal- 
licity and age would have t o be coordinat ed such that 
(AZ)/Z = -0.7(AAge)/Age (| WorthevI [1993 ) . In fact, an 
anti-correlation with exactly such a slope is indeed observed 
from spectr oscopic analyses which break the age-m etallicity 
degeneracy (|Trager e t al."2000l: [sinith et al. ' 2008a'). So it is 
possible that there is greater variation in the stellar popu- 
lations of bulges than their colours suggest. 



the error bars on the slopes are considerably (a factor 2-3 
times) larger and so no firm conclusions can be drawn. 



6.2 Density as an Environment Indicator 

The cluster-centric radius is a reasonable proxy on average 
for the depth of the cluster potential well , but in reality clus- 
ters are built from groups and singles (see lMcGee et al.ll2009l . 
for a detailed discussion.). If the environmental influences 
are attributable to a given galaxy's previous ("group-like") 
environment ( "pre-processing" ) then an environmental indi- 
cator based on the local density might be more sensitive to 
such effects. 

To test this we define the local density E5 as follows. 
Using galaxies brighter than our absolute magnitude limit 
of —19.5, we find the distance to the fifth nearest neigh- 
bor in the cluster and compute the density using this ra- 
dius. It turns out that this definition is almost identical 
(w ithin a few hundredt hs of a magnitude) to the limit used 
by iBalogh et al.l (|2004l '). who used SDSS r mags. However, 
rather than defining the density within a ±1000 km s~^ red- 
shift slice, we simply consider only galaxies that are clusters 
members. 

As seen in Figure (TT] the total colour varies as 0.032 ± 
0.008 per dex in density, in the expected sense that galaxies 
in higher density environments are redder. The bulge colour 
slope is not statistically significant (0.006 ± 0.019) but the 
disc colour slope is 0.050 ± 0.013, significant at the 3.9 a 
level. The significance of the trends with density are similar 
to, or slightly lower in significance than the radial trends as a 
function of log(rp/r2oo), suggesting that these two measures 
of environment are largely degenerate when applied to rich 
clusters. 



7.2 Disc fading 

The dependence of disc colour on radius strongly sug- 
gests a scenario in which their star formation has been 
shut off by the cluster environment, after which the discs 
fade and redden. It is interesting to determine whether 
this process is abrupt, as one might expect if the cold gas 
is stripped by ram pressure, or slow er, as might be ex - 
pected in the "strangulation" scenario (jLarson et al.lll980l ). 
in which the hot gas halo is stripped, but the cold gas 
remains to f uel star formation . The former scenario cer- 
tainly exists: ISmith et al.l (|2010h found that at least 1/3 of 
blue galaxies in the 1 /i~^Mpc-radius core of Coma show 
UV tails extending from the disk suggestive of star for- 
mation during ram pressure stripping, and the direction of 
the tails suggests that the stripping is occurring on infall. 
It is not clear, however, how efficient and how rapid this 
process is. Evidence for a slower "strangulation" mecha- 
nism is favour ed by a co mpariso n of models with red and 
blue fractions (JFont et al. 2008: van den Bosch et al.ll2008l : 
IWeinmann et al.ll2009l : iMcGee et allbOOgH . 

We have found that the difference in median colour be- 
tween discs at the centre of the cluster and those at the virial 
radius is 0.103±0.023 magnitudes in B — R. It is interesting 
to see if this colour difference can be reproduced with simple 
infall models in which star formation is truncated as the disc 
falls into the cluster. We use the simulations of ICao et al.l 
(2004) to associate a given cluster-centric radius with a me- 
dian infall time. Their Figure 15 suggests that while objects 
in the centre of the cluster were accreted at 2; = 1 (a look- 
back time of 7.7 Gyr), those at the virial radius typically 
were accreted at z ~ 0.4 (4.3 Gyr ago). 

If we assume that solar metallicity discs began form- 
ing stars at 10 Gyr ago at a continuous rate until cross- 
ing the cluster virial radius, followed by abrupt quenching, 
we obtain a colour of 1.43 for the disc in the centre of the 
cluster, in reasonable agreement with the median observed 
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Figure 17. Total colour, bulge colour, and disc colour versus projected density. 



value of 1.38. However, at the virial radius, the predicted 
colour is 1.40, which is quite red compared to the observed 
colour of 1.28. Put another way, the predicted colour differ- 
ence (a more robust quantity than the absolute colour) is 
0.03 whereas the observed value is 0.10 ± 0.02. If we instead 
adopt a "strangulation" scenario for the disc in which star 
formation is not cut off abruptly, but instead declines with 
an exponential decay rate of 1 Gyr after infall, the predicted 
colour difference between a disc at the centre and one at the 
virial radius is 0.07. This is in better agreement with the 
observed radial gradient in disc colour. 

However, the above naive interpretation is significantly 
complicated by projection effects. There are two factors at 
play: first, an intrinsic three-dimensional colour gradient will 
be flattened when viewed in projection; second, there will 
be contamination by infalling field galaxies that are outside 
the virial radius but which appear within the virial radius 
in projection (and which have sufficiently low velocities to 
pass the ±3aci clipping). These field interlopers are found 
in the outer regions and so will steepen the observed colour 
gradi ent. The fraction of su ch interlopers is difficult to deter- 
mine. |Mamon_et_alJ (120101 ) estimate that the fraction of such 
interlopers is a strong function of projected radius. More 
specifically, for the subsample with r-p > 0.3r2oo, one might 
expect ~ 33% of objects to lie beyond the virial radius in 
three dimensions (G. Mamon, private communication). As 
a further complication, most of these "interlopers" will be 
infalling blue galaxies, but some may be "backsplash" galax- 
ies which have passed through the centre of the cluster and 
have exited the virial radius. 

We can test the null hypothesis that there is no intrin- 
sic variation within the cluster virial radius. Specifically, we 
simulate the colour distribution of the outer radial bin with 
a mixture of the disc colours from the "core" bin (67%) with 
33% contamination from the field discs. While the median 



colour of this simulated population is slightly redder than 
the observations, a K-S test shows that there is no signifi- 
cant difference between it and the observed distribution of 
disc colours in the outer radial bin. Thus we cannot rule out 
the hypothesis that the reddening begins beyond the virial 
radius. Nevertheless, the data do not exclude some redden- 
ing within the cluster itself: indeed a better fit is obtained 
if we assume that, intrinsically, the core galaxies are red- 
der than the outskirt galaxies by 0.03 mag (in rough agree- 
ment with the "abrupt quenching" model discussed above). 
In summary, given the uncertainties associated with the in- 
terloper contamination fractions, and with the small size of 
the comparison field sample, it is difficult to decontaminate 
the cluster samples with sufficient accuracy in order to de- 
termine which scenario (abrupt quenching or strangulation) 
is at play. Nevertheless, we emphasize that there is indeed 
a difference between the disc colours of cluster-core galaxies 
and those at larger radii. 



7.3 Implications for the Morphology-Density 
Relation 

In the previous section, we showed that the data are con- 
sistent with a scenario in which the disc component fades 
within the cluster environment. Is this disc fading consis- 
tent with the buildup of the morphology-density relation? 
More specifically, are SOs formed by spirals after their discs 
fade and their spiral density waves are damped? 

To address this question, we first need to estimate the 
amount of fading. Comparing a hypothetical field spiral 
where the disc continues forming stars at a const ant rate to 
one w hich was truncated 6 Gyr ago, the models of lMarastonI 
((2OO5I) predict a difference of a factor of ~ 2.5 in the 7?-band 
stellar mass-to-light ratio and a factor ~ 4 in the B-band. 
Thus, assuming the bulge is unaffected by environment, as 
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a galaxy's disc fades, the total R luminosity will drop by a 
factor less than 2.5, and B/T will increase by a factor up 
to 2 .5. General support for this picture conies from stud- 
ies IjNeistein et al. 1993: Bedregal et alj 12006) which have 
shown that the So ItuUv fc Fished l|l977l . hereafter TF) rela- 
tion is approximately consistent with being a faded version 
of the spiral TF relation. 

However, there are problems with applying this simple 
fading scenario to map the field population into the cluster 
population as a whole. One problem arises when comparing 
the cluster B/T distribution to that in the field. In Fig.[8l we 
showe d the field B-band B/T distribution from I Allen et al.l 
l|2006ri . which has a large excess of late types with B/T < 
0.1 (and many with no bulge at all). As noted above, in 
the _B-band, the fading might be as much as a factor of 4 
between the field and the cluster. Thus field late- type spirals 
which would have had B/T < 0.1 had their star formation 
continued, might appear as B/T ~ 0.2 if their star formation 
was truncated at cluster infall. However, Fig. [8] does not 
show an excess of objects with B/T ~ 0.2, but rather a 
broad distribution peaking at B/T ~ 0.5 in the B-band. 
Thus in order to create the cluster B/T histogram from 
that in the field, one needs to either transform the discs 
of late-type spirals into bulges, perhaps by "harassment" 
l|Moore et al.lll998h . or to remove the very low B/T objects 
from the distribution entirely, either by fading them below 
the magnitude limit or by destroying these galaxies. 

A further problem with the scenario in which all SOs are 
formed from present-day spirals is the observation that the 
bulges of SOs are larger and brighter than those of spirals 
llDressleilll980l : iBurstein et"aLll2005l : IChristlein fc Zabludofj 
[2OOJ). Therefore it seems that while the brighter SOs in 
present day clusters are unlikely to be the faded remnants 
of today's disc galaxies, this does not exclude fading as pos- 
sible mechanism for creating lower-lumtnosity SOs. Whether 
the brighter cluster SOs are consistent with being the faded 
versions of higher-redshift spirals remains to be seen. 



8 CONCLUSIONS 

We have shown that the "tilt" of the total-light CMR is 
driven primarily by the tilt of CMR of the bulges (but not 
by the discs), and to a lesser extent, by a change in the B/T 
ratio as a function of magnitude along the red-sequence. 
The tilt of the CMR is therefore due primarily to stellar 
population changes in the bulge as a function of magnitude 
or a. 

In clusters, we find that discs in L, galaxies are ~ 0.25 
magnitudes bluer than bulges. If we assume that discs have 
solar metallicity, their colours are consistent with an SSP 
age of 4.5 Gyr. 

Although bulge colours are insensitive to environment, 
the disc colours show a statistically significant dependence 
on cluster-centric radius in the sense that discs are bluer at 
the virial radius by ~ 0.103 ± 0.023 magnitudes in B ~ R 
compared to their counterparts of the same total magni- 
tude in the cluster centre. This effect drives the total-colour- 
radius relation noted by previous authors. A straightforward 
interpretation of this result is that it is due to the quenching 
of star formation in the disc upon cluster infall. 

In summary, the physical processes that affect star for- 



mation in bulges are mass- dependent, and presumably pri- 
marily internal, whereas star formation in the disc shows 
little dependence on mass, but a strong dependence on en- 
vironment as parametrized by cluster-centric radius. 
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APPENDIX A: AGES AND METALLICITIES 
OF BULGES AND DISCS FROM COLOURS 
AND SPECTRA 

Bulges and discs have different colours, presumably due to 
differing stellar populations. Due to degeneracies, it is dif- 
ficult to measure both ages and metallicities from optical 
colours alone. Spectra allow one to break this degeneracy, 
and ideally, one would use IFU spectra to obtain spatially 
resolved ages and metallicties, but, at present, such IFU 
samples remain small and somewhat het erogeneous (e.g. 
iRawle et al] l2008l : iKuntschner et all [2O1O). However, the 
NFPS contains spectra and stellar absorption linestrengths 
for 4097 red galaxies obtained from 2 arcsecond diame- 
ter fibres. From the linestrengths, a strong "downsizing" 
trend - younger a ges in less massive galaxies - was found 
l|Nelan et al.ll2005l ). However, if the disc is younger than the 
bulge, and if B/T varies along the red-sequence, some of 



this trend may be due to variations in the fraction of central 
spectra arising from discs. 

It is possible to use these spectra to infer bulge ages if we 
make some assumptions about the disc properties. Since the 
disc colour is independent of magnitude and a, we will as- 
sume that discs have an SSP age of 4.5 Gyr, solar metallicty 
and solar ratios of g-eleme nts to iron. These parameters are 
chosen becau se they yield Maraston! (jl_998) model colours 



(assuming a (IKroupaluOOl ) initial mass function), for the 
disc consistent with their median observed colour of 1.353. 
We will then model the contributions to the spectra from 
the bulge and disc components, and fit the stellar popula- 
tion parameters of the bulge using the central spectra. We 
will do this for subsamples binned by velocity dispersion, 
rather than magnitude because it has been shown that ve- 
locity disper sion, a, is th e driving parameter of the stellar 
populations (jSmith et al.. 2009a; Graves ct al. 2 009 ) . 

The central 2 arcsecond aperture contains a mixture of 
bulge and disc light, and so careful modelling is needed. In 
order to assess the contribution from bulge and disc to the 
central spectra, we use the GIM2D parameters to first calcu- 
late the bulge and disc light within a 2 arcsecond diameter 
aperture (matching the NFPS spectra), assuming that the 
exponential disc light can be extrapolated to the center. We 
have found that, not only is a strongly correlated with stellar 
populations, it is also strongly correlated with morphology: 
Fig. lAll shows B/T and the aperture B/Tap, as a function 
of a. Note the strong correlation and relatively low scatter 
(compare with the weak trend and high scatter in B/T ver- 
sus magnitude in Fig. |6}. Nevertheless, at all \og{a) > 1.8, 
the central light is always dominated by the bulge. Thus for 
each a bin, we know B/Tap from Fig. lAll 

We proceed to model the ages and metallicities as fol- 
lows. Given the disc age and metallicities, we vary the bulge 
age, metallicity and a- enhancement. For each choice of bulge 
parameters, we calculate the stellar mass-to-light ratios of 
the bulge and disc components and, given the fraction of 
light in the bulge from B/T^ , we solve for the stellar mass 
fraction in the bulge and in the disc comp onents. We use 
the co mposite stellar population mo dels of Allanson et al.l 
1120091) ■ which in turn are based on iThomas et al.l (120031 . 
120041 ). to calculate spectral linestrength indices for the two- 
component model and compare these to the mean NFPS 
linestrengths for the same a bin, and then minimize to find 
the best-fit bulge age and metallicity. 

The resulting bulge ages and metallicities as a func- 
tion of a are shown in Fig. IA2I We show results for two 
data samples: the NFPS (Nel an et al.|,2005|); an d 3 clusters 
in the Shapley concentration (|Smith et al.ll2007l ). which are 
also observed through a 2 arcsecond fibre, and are at a sim- 
ilar distance to the NFPS clusters studied in this paper. In 
both cases, the effect of including a 4.5 Gyr disc compo- 
nent is to shift the bulges to older ages, typically by 1 — 2 
Gyr, compared to simple single-component SSP fits. There 
is some disagreement between the two samples for the low- 
est logo- ~ 1.8 bin, for which the NFPS requires bulge ages 
~ 4.5 Gyr, whereas the Shapley sample prefers ages ~ 6.5 
Gyr. For both samples, the "downsizing trend", in which 
bulges with larger velocity dispersions are older, is present 
in bulge components. 

Of course, this result is based on the assumption that 
disc surface brightness profile continues to be an exponen- 
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tial right to the centre of the galaxy. If the surface bright- 
ness profile flattens towards the centre, as is common in 
SO's (|Freemanlll97d ). then the disc contribution is reduced 
compared to the bulge, and our conclusion regarding bulge 
downsizing is strengthened. It is, of course, in principle also 
possible that disc age does vary along the sequence, and 
that metallicity and/or dust conspire to achieve a disc colour 
which is independent of magnitude or a, but such scenar- 
ios would require fine-tuning and so seem rather contrived. 
Thus we conclude that the variation in bulge metallicity 
along the red-sequence appears to be relatively less impor- 
tant than bulge age in driving the bulge CSR (and hence 
CMR). 
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Figure Al. ij— band B/T as a function of velocity dispersion, cr. The left panel shows B/T, whereas the right panel shows the B/Tap, 
the B/Tobtained from the light within a 2 arcsecond diameter aperture, chosen to match fibre spectroscopy. As in other plots, the solid 
squares show medians fo bins of a. Note the tight correlation between log a and B/T. 





Figure A2. Bulge age as a function of bulge metallicity. The left panel shows fits to the NFPS central linestrength data, assuming a 
4.5 Gyr solar metallicity disc (indicated by the solid circle). Contours indicate 1 and 2 a confidence levels on two parameters jointly 
for bins with (log cr) = 1.94, 2.09, 2.17, 2.26, 2.39 from b o ttom left to top right. The base of each arrow shows the SSP parameters if we 
assume only a single component SSP as in lNelan et al.l 1120051) . whereas the tip of the arrow shows the best fit bul ge parameter s of th e 
bulge+disc model. The right panel is the same except for the data fr om 3 clusters in the Shapley Concentration from lSmith et al.l 1(2003), 
in comparison to single-component SSPs from I AUanson et al.l 1120091 ) for bins (logcr) = 1.71, 1.84, 2.02, 2.16, 2.31, from leftmost to upper 
right. Note that the central spectra require that the bulge component of lower a galaxies is younger than that found in larger a galaxies 
(bulge "downsizing" ) . 



